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The halo factor is one of the experimental data which describes a distribution of neu- 
trons in nuclear periphery. In the presented paper we use Skyrme-Hartree (SH) and 
the Relativistic Mean Field (RMF) models and we calculate the neutron excess factor 
A b defined in the paper which differs slightly from halo factor / CX p • The results of the 
calculations are compared to the measured data. 

1. Introduction. 

One of the important problems of nuclear physics is the description of the periphery 
of the nucleus. From theoretical and experimental point of view the topic is com- 
plex and much information is needed in order to predict neutron or proton excess 
in the neighbourhood of nuclear surface. 1 ' 2 ' 3 The peripheral density distributions 
of protons and neutrons are determined theoretically by single particle eigenfunc- 
tions of protons and neutrons in this region. Applying different nuclear models one 
calculates wave functions in the deformed oscillator basis. They behave asymptoti- 
cally like exp(— 7r 2 ), where 7 is a constant and r is the distance from the center of 
the nucleus and do not reproduce properly the particle densities at large distances. 
The real nuclear densities behave like exp(— (3r). In the following we rather stay in 
the nuclear interior where the density distributions are well defined and we describe 
the nuclear peripheral properties using only well defined theoretical internal nuclear 
densities. We would like to show that a simple function 8b {r) introduced by Bethe 5 

&b{t) = j^PN(r) - pz(r) = Pn{t) - pzif) , (1) 

allows to predict the experimentally observed halo data without the knowledge of 
the local behaviour of a tails of nuclear densities in asymptotic region. The function 
given by Eq.(l) is the difference of the relative neutron density pn{t) = (Z/N)pn 
and the proton density pz(f) at a point r. 

In the following section we introduce a neutron excess factor (NEF) based on 
5b defined in Eq. (1) function and we calculate NEF for a class of nuclei studied at 
Low Energy Antiproton Ring (LEAR) facility in CERN. 1,2 . Two mean field models 
are used in calculating NEF which after that is compared to the experimental data. 
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Fig. 1. Bcthe function of 208 Pb (multiplied by r 2 , solid line) and the integral Ib(x) (dashed line). 
The maximum of Ig corresponding to R = 6.5fm (zero of 5g(r)) is the neutron excess factor 
discussed in the text. 



2. Model 

It is a trivial to notice that the integral of the Bethe function (1) over the whole 
space vanishes 

J S B {r)d 3 f=0. (2) 

This direct property allows to define the total peripheral neutron or proton excess 
without the knowledge of the density distribution in the nuclear periphery. For 
example, in the case of the spherical symmetry equation (2) reads 

/>oo pr' f-oo 

4ir r 2 5 B (r)dr = 4tt r 2 S B (r) dr + 4tt r 2 6 B (r)dr = 0, (3) 

JO JO Jr> 

where r' is an arbitrary point which if chosen properly will maximize the first term 
in the above equation. Let us define the following function of the distance r 

I B (r) = -4tt / r' 2 dr' . (4) 
Jo 

The maximal value of A B with respect to r 

A B = m&x r I B (r) , (5) 

gives the relative outer excess of neutrons or protons out of the point which max- 
imizes I B . The positive value of indicates the excess of neutrons while the 
negative value of it corresponds to the case of the relative proton excess. In the 
following we call the value of the neutron excess factor (NEF). Both S B and B 
functions for 208 Pb are shown in Figure 1. 

Similar functions and quantities are defined for more general case of a deformed 
nuclei. 4 The corresponding integrals and quantities (Eqs. 4 and 5) are calculated in 
this case by taking the integrals over the region of the deformed surface similar to 
the average nuclear surface. 
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Fig. 2. Annihilation of pp or pn on the periphery of the deformed nucleus. The inner and the 
outer regions are divided by a sphere of radius r. The pions are the outcome of the annihilation 
reaction. 

3. Antiproton annihilation on nucleus 

The measurements of products of the annihilation of antiprotons on nuclei 1,2 inform 
on the value of the neutron excess in the nuclear periphery where the annihilation 
takes place. This is the outer region extending over the whole periphery of the 
nucleus. From the experiment 1,2 one approximatelly knows the average radius of 
stopping of antiprotons. It is slightly larger (by about 2fm) than the average nuclear 
radius R = r yl 1/3 . 

While the deformed nuclei are not oriented in the target the information about 
the process is contained in the averages which depend neither on the direction 
nor on the distance from the center of the nucleus (see Fig. 2). This corresponds 
to the averaging over the outer region. Assuming the annihilation cross section 
locally proportional to the density of nucleons (protons or neutrons) one sees that 
the ratio of the cross section of stopping the antiproton on neutrons to the cross 
section of stopping antiproton on protons is proportional to corresponding ratio of 
the densities pn/pz- The total annihilation ratio is thus approximately proportional 
to the integral of the ratio of the densities over the whole outer region. A similar 
averaging has been performed in calculations of and therefore the result given 
in Eq. (5) i.e., NEF is proportional to the halo factor / oxp (see Eq. (6)) defined 
by 1,2 

f „ — ^ (Q) 

/cxp N Afp p ' W 

where Np p and Np n denote the number of annihilations of both pp and pn types 
respectively (the final number of Z -\X N nuclei corresponding to pp annihilation 
and the final number of zX" -1 nuclei which correspond to pn annihilation of p on 
a nucleus zX N ) The final nuclei are "counted" after the anihilation processes has 
taken place. 

4. Results 

In the following section we compare our results of calculations of As to the exper- 
imentally measured halo factor / oxp . 
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Fig. 3. Neutron excess factor Ag (NEF) for nuclei studied by the Low Energy Antiprotonic Ring 
(LEAR) group at CERN. Experimental data are given with crosses connected by solid line. Plain 
circles represent the RMF data and the open circles (dashed line) correspond to Skyrme-Hartrcc 
calculations with SkVII force. For more details and the references see the text. 



Figure 3 shows the outer neutron excess factor (As) as calculated in the Rel- 
ativistic Mean Field model (RMF, plain circles) with NL3 set of parameters 6 and 
in the selfconsistent Skyrme-Hartree (SH) model with SkVII force (open circles). 7 
The data are presented for nuclei reported by Low Energy Antiprotonic Ring group 
(LEAR) at CERN. 1 ' 2 It is observed that the calculations in both RMF and SH mod- 
els arc strongly correlated with experimental halo factor / oxp (the scale of As is 
chosen arbitrarilly) . The correlation shows that the NEF Ab is a good indicator of 
the neutron halo property and it gives the measure of the peripheral annihilation 
ratio of the antiproton on the neutrons and protons in the nuclear system. 
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